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ABSTRACT 

We consider current observational constraints on the presence of cool, optically thick disk material in 
quiescent black hole binaries, specifically focusing on a case study of the prototypical system A 0620- 
00. Such material might be expected to be present theoretically, but is usually claimed to make 
a negligible contribution at optical and infrared wavelengths. The primary argument is based on 
measurements of the veiling of stellar photospheric absorption lines, in which it is assumed that the disk 
spectrum is featureless. We use simulated spectra to explore the sensitivity of veiling measurements 
to uncertainties in companion temperature, gravity, and metallicity. We find that the derived veiling 
is extremely sensitive to a mismatch between the temperature and metallicity of the companion and 
template, but that the effect of a plausible gravity mismatch is much smaller. In general the resulting 
uncertainty in the amount of veiling is likely to be much larger than the usually quoted statistical 
uncertainty. We also simulate spectra in which the disk has an emergent spectrum similar to the 
star and find that in this case, optical veiling constraints are moderately robust. This is because 
the rotational broadening of the disk is so large that the two line profiles effectively decouple and 
the measurement of the depth of stellar lines is largely unbiased by the disk component. We note, 
however, that this is only true at intermediate resolutions or higher, and that significant bias might still 
affect low resolution IR observations. Assuming that the optical veiling is reliable, we then examine 
the constraints upon the temperature and covering factor of any optically thick disk component. 
These are stringent if the disk is warm (T e g > 3500 K), but very temperature sensitive, and cooler 
disks are largely unconstrained by optical measurements. Current IR veiling estimates do not help 
much, representing rather high upper limits. Probably the best constraint comes from the relative 
amplitudes of ellipsoidal variations in different bands as these are sensitive to differences in veiling 
which are expected for disks cooler than the companion star. A significant disk contribution in the 
IR, up to ~ 25 %, is not ruled out in this or any other way considered, however. 
Subject headings: accretion, accretion disks — binaries: close 



1. INTRODUCTION 

Almost all of our knowledge of the masses of stellar- 
mass black holes derives from those in low mass X-ray 
binaries: the black ho le X-ray transients (BHXRTs; e.g. 
iCherenashchuld 1200(1) . In most of these systems sub- 
stantial X-ray activity occurs only during well defined 
outbursts. Between outbursts, the emission from the ac- 
cretion flow fades to the point that the companion star 
is clearly visible and nearly undisturbed by irradiation; 
hence it can be used to derive a dynamical mass for the 
black hole. Current methodology relies on three key in- 
gredients to determine a mass. Firstly the radial velocity 
curve of the companion star yields the binary period (if 
this is not already known), and the companion's pro- 
jected radial velocity semi-amplitude. Based on these 
quantities alone, the mass function can be calculated, 
and this sets an absolute minimum mass for the compact 
object. Secondly the rotational broadening of the com- 
panion star's absorption lines is used to derive the mass 
ratio. Finally, using photometric observations (prefer- 
ably in the IR) of the ellipsoidal modulation of light from 
the companion star, coupled with models of this modula- 
tion the binary inclination can be determined as a func- 
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tion of assumed mass ratio. Together these quantities 
yield the mass of the black hole. 

In the ideal case there will be no light at all from the 
accretion flow. This is never completely the case, how- 
ever, and we usually must make do with making some 
assumptions about its effect on the analysis. Typically it 
is assumed that the spectrum is a featureless continuum 
and that its contribution is minimized in the IR. These 
assumptions are q uestionable, h owever. The disk insta- 
bility model (DIM: lLasotal200ll) advanced to explain the 
transient behavior predicts that the quiescent disk should 
be very cool, likely cool er than the companion star, and 
could be optically thick l)Cannizzoll993tlMenoul200"2T) . It 
has been argued that such a disk is virtually required to 
account for the necessary surface densities, although we 
argue in Section that this requirement is much weaker 
than usually stated. Under these conditions we might 
expect both a significant IR contribution and possibly 
absorption lines similar to those of a late-type star. Un- 
accounted for disk emission will bias inclination deter- 
minations in the sense that inclinations will be under- 
estimated if disk light is not taken into account. This 
effect is quantified by measuring the dilution of the pho- 
tospheric absorption lines, and is almost always done in 
the optical. One must then ask how robust this estimate 
of the dilution (the veiling) is, whether it is biased by 
absorption lines from the disk, and whether the veiling 
from a cool disk could be significantly larger in the IR 
than in the optical. 
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We know that a disk is present in quiescence. This 
is indicated by line profiles and Doppler tomography 
of emission lines (e.g. iMarsh. Robinson, fe Woodl I1994J) 
and probably by observation s of flic kering behavior (e.g . 
iZurita. Casares. fc Shahbal I2003T: iHvnes et alJ 12003ft . 
The line observations, at least, are related to hotter 
and optically thin gas which probably co-exists with the 
cool optically thick material that we are considering, ei- 
ther in a chromospheric layer above the disk, or in low 
surfa ce density regions ijVrielmann. Hessman. fc Hornd 
2002). The origin of the flickering is less clear. In 
at least one case, V404 Cyg, it has a line component 
which appears to be d istributed across the d isk and 
driven by X-ray heating l)Hvnes et alJ2002ll200^1 . These 
flares exhibit a continuum as well, but this too may be 
optically thin, and/or associated with the inner disk. 
None of these observations inform us about a possible 
cooler optically thick component of the disk. We also 
note that in the quiescent white d warf system, WZ Sge, 
IHowell. Har rison. fc Szkodvl {2004) have recently identi- 
fied CO and H2 in emission. The authors argue that this 
originates in a cool (2000-5000 K) component of the disk. 
No equivalent detection has yet been made in quiescent 
BHXRTs, and WZ Sge is a unique object that may not be 
comparable to them, but this should be borne in mind. 

In this work, we will consider what other evidence ex- 
ists to constrain such cool material. We will focus on 
one system as a case-study, the prototypical BHXRT, 
A 0620-00; this work is intended to be an example of the 
severity of different sources of uncertainty rather than an 
exhaustive exploration of parameter space. We will begin 
by reviewing the relevant system parameters. We then 
construct synthetic line profiles, and model spectra, for 
both the companion star and an optically thick disk with 
a cool, stellar spectrum. Then we use these to test how 
sensitive veiling measurements are to uncertainties in the 
temperature, gravity, and metallicity of the companion, 
and whether the veiling measurements may be biased 
by the presence of such disk absorption lines. Assuming 
that the optical veiling estimates are reliable, we next 
examine the constraints this places upon the quiescent 
disk. Finally we will examine potentially more stringent 
constraints from the IR veiling and photometric color. 

2. A 0620-00: PARAMETERS FOR A CASE STUDY 

A 0620-00 i s one of the best studied syste ms with 
both optical llMarsh. Robinson, fc Woodl Il994|) and IR 
(Shahbaz. Bandvopadh vav. fc Charlesl 119991) veiling es- 
timates a s well as recent high-resolution optical spec- 
troscopy ((Gonzalez Hernandez et al1l2004[) . 

The temperature of the companion star is criti- 
cal for veiling estimates, but is typically very uncer- 
tain in BHXRTs. Several of the methods used may 
be compromised if the disk contamination is signif- 
icant, and does not vary monotonically with wave- 
lengt h. In particular, photometric sp ectral types (e.g. 
K4V, iGelino. Harrison, fc Oroszll200l|) could be signifi- 
cantly biased by contamination of the continuum b y an- 
other source of light, a point also appreciated bv lOkei 
l|1977|) . Si milarly spectr ophotometric constraints (e.g. 
K3V-K4V; Haswcll. 119921) have the same problem. Con- 
straints based purely upon spectral features are rela- 
tively weak for A0620-00, with the excep t ion of the re- 
cent work by IGonzalez Hernandez et all <)2004[) which 



we will discuss later. Mos t have use d broad bands 
rather than sharp features. lOkei l)1977fl obtained K5- 
K7, and preferred a dwarf to a giant. Subsequent obser- 
vations have found molecular bands to be weaker, how- 
ever, and allow earlier types. iMcClintock fc Rlem illard 
(1986) used K2V, K3V, K4V, and K7V t emplates 
and d id not explicitly reject any of them. lHaswelll 
(1992) preferred K3V or K4V, but could not re- 
ject K8V. Earlier spectral types are possible as the 
only other type considered (an d rejected) was G8V. 
IMarsh. Robinson, fc Wool < 1994ft only used a K3V-K4V 
templ ate" Finally IsTiahbaz. Bandvopadhvav. fc Charlesl 
( 1999) examined the CO bandhcad in the K band and 
compared to K0V, K3V, K5V, and K7V spectral types. 
K3V was preferr ed, but nearby types could n ot be ruled 
out, especially as lFroning fc Robinson! l|2001f l challenged 
the validity of results based purely on the CO bandhead 
because of its sensitivity to surface gravity and metal- 
licity. In summary, then, most studies allow a range of 
spectral types. A classification as early as K2V, or per- 
haps even K1V, cannot be definitely rejected. Similarly 
later types, K4V-K5V, perhaps as high as K7V-K8V are 
not implausible considering only spectroscopic evidence. 
For this work w e will assume K4V, and a temperature 
of T cff = 4690 K lIGrav fc Corballvlll994ft . We have bela- 
bored this point more than strictly necessary to illustrate 
the uncertainties present eve n for a well studied object. 

The more recent analysis bv lGonzalez Hernandez et alJ 
l)2004ft is superior to those discussed above, and yields a 
somewhat hotter temperature, T c g — 4900 ± 100 K, cor- 
responding to a spectral type of K3V. We deliberately do 
not use this value in our simulations as it is not repre- 
sentative of the quality of measurement available in other 
systems. We do keep it in mind as a further indication 
that typical temperature estimates are uncertain by at 
least 200 K (and often more) and will comment further 
where appropriate. 

We also require a surface gravity. We a dopt the 
predicted value of IGelino. Harrison, fc Pros! 1)200 lj) of 
log g = 4.46 ± 0.04, based on the derived system pa- 
rameters. It is consistent with the measured (spec- 
troscopic) value of 4.2 ± 0.3. The small uncertainty 
in log q is dominated by t hat i n the mass ratio of 
IMarsh. Robinson, fc Wood lfl99l of q = 0.067 ± 
0.010. If the inclination is large r than that derived by 
IGelino. Harrison, fc Orosa (|200 IT) , for example if IR veil- 
ing is significant, then the effect on log g is relatively 
small. For i = 50° and 60° it decreases to 4.40 and 4.35 
respectively. 

Finally for the companion we require a metallicity 
and rotational broadening velocity. For the metallicity 
we assume solar, since this is the default assump tion 
for most systems. IGonzalez Hernandez et aTl l|2004ft ac- 
tually find slightly super-solar metallicity (0.14 ± 0.2), 
but this is consistent with solar within uncertainties. 
For the rotational broaden i ng we adopt 83kms _1 from 
IMarsh. Robinson, fc Woodl l|1994ft . 

To compute disk line profiles we additionally require 
orbital parameters. For this purpose we assume P or b = 
7.78 hrs, q = 0.067, M, = 11 M^, and i = 40.75° 
( Geli no. Harrison, fc Oroszl2001l and references therein) . 

Veiling measurements hav e been made in both opti- 
cal an d IR for A0620-00. IMarsh. Robinson, fc Woodl 
((1994ft measured disk contributions of 6 ± 3% at 
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Fig. 1. — Examples of composite line profiles calculated in the companion star's rest frame for various resolutions. We have used 
parameters appropriate to A0620— 00, except that we have increased the disk contribution to 50%, so that star and disk contribute equally. 
The solid line shows the profile derived if the disk has a line of identical equivalent width to that of the companion. The dashed line 
assumes the disk is a featureless continuum. Resolutions are chosen to indicate cases cor responding to an arbitrarily high resolution 
(R = 30 000), a typical optical resolution (Ft = 4000, e.g. IMarsh. Robinson. Kr. Woodlll994fl . and a typical IR resolution (R = 500, e.g. 
Shahbaz. Bandvopadhvav. & Charles 1999). 



Ha an d 17 ± 3% at Hff usin g a K 3-K4V tem- 
plate. IGonzalez Hernandez et aTl l|2004|) found com- 
parable values (0 ± 5% and 20 ± 5% r espectively). 
iShahbaz. Bandvopadhvav. fe Charles! (|1909f) obtained an 
upper limit on the disk contribution in the IR of 27% 
(2 — a). This is not a robust result, but it is all we cur- 
rently have. Over and above the limitations discussed 
explicitely later in this work, two major concerns are i) 
the data quality is not high enough to be confident that 
the main CO features are securely detected, and ii) cat- 
aclysmic variable star secondaries exhibit anomalously 
weak CO features, po ssibly due to prior CNO processing 
l|Harrison et alJl200l . The same mechanism may oper - 
ate in BHXRT companion stars (?; lHaswell et~atl [2002). 
calling into question the usefulness of CO bands for veil- 
ing measurements. 

3. THE SIMULATIONS 

3.1. Methodology 

In general, our analysis technique is to create model 
composite spectra including a rotationally broadened 
model atmosphere spectrum for the companion, and ei- 
ther a flat continuum or a model atmosphere for the disk, 
convolved with a disk line profile. Throughout we aim 
to reproduce the main characteristics in order to obtain 
illustrative results; we are not attempting precise, state- 
of-the-art modeling of the quiescent spectra. The syn- 
thetic spectrum can be expressed as (l — a)F atax +aFdi$k, 
where a is the assumed true disk fraction (veiling) , F s t a r 
is the input stellar spectrum, and Fdi s k ^ s either unity 
everywhere, or an input disk spectrum. Having con- 
structed the model spectrum in this way, we then fit it 
with (1 — /3)i 7 tcmpiatc + where (3 represents the derived 
veiling and -^template is the spectrum of the adopted tem- 



plate star. Note that in fitting templates to model data 
we always assume the disk is a flat continuum, following 
observational practice. The difference between a and (3 
then indicates how much the veiling is biased by choice of 
template, or the presence of disk spectral features. Our 
treatment is intended to represent the ideal observational 
case, so we use noise-free normalized spectra. In ana- 
lyzing actual observations further uncertainties will be 
introduced by noise and difficulties in normalizing the 
spectra. 

We use two representative wavelength ranges. For the 
optical we use 5920-6520 A; this excludes the strong NaD 
lines as these might be contaminated by interstellar ab- 
sorption in observational data, and also avoids the Ha 
and He I 5875 A emission lines. In the IR we use 2.19- 
2.31 /mi, to exclud e Br7, but include the first CO band- 
head (as used by IShahbaz. Bandvopadhvav. fc Charles! 
1999). We calculate models appropriate to three in- 
strumental resolutions, represented by Gaussians. The 
reference case uses R — A/AA = 30 000, sufficient 
to resolve the stellar rotational broadening profile. 
The second case is R = 4 000, representative of op- 
tical studies, and based on th e specific example of 
IMarsh. Robinson, fc Wood! l) 19941) . The third case is R = 
50 0, based on the IR veiling determination for A 0620-00 
bv IShahbaz. Bandvopadhvav. fc CharTesI <)1999|) . 

3.2. Synthetic Line Profiles 

For convenience we calculate all model line profiles in 
the rest frame of the companion star; in practice ob- 
servers often shift spectra into this frame and average, 
so this is appropriate. 

We represent the stellar line profile with a rotational 
broadening profile ijGravl Il992;) assuming linear limb- 
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darkening with a limb-darkening coefficient of 0.73 in the 
opt ical, appropriate around 6000-6500 A and 0.31 in the 
IR |)A1-Naimivlll978() . This is obviously a crude approx- 
imation to a distorted Roche-lobe filling star, but is ad- 
equate for our purposes. Likewise, linear limb-darkening 
is imprecise, but commensurate with other approxima- 
tions made. 

Disk profiles assume a disk extending from i?j n to i? ut- 
i?i n is taken to be 10 4 Schwarzschild radii, but in prac- 
tice has negligible effect on the line profiles. We expect 
i? ou t to lie between the circularization radius, the mini- 
mum radius at which new material is likely to be added 
to the disk, and the tidal truncation radius, the approxi- 
mate maximum size. The former may be more likely for 
a quiescent disk, but the latter produces narrower pro- 
files, and hence is more likely to contaminate the stel- 
lar spectrum. For A 0620-00, these assumptions yield 
v ou t smi = 600kms _1 for the tidally truncated case and 
v ou t sin i = 740 km s _1 for a disk extending out to the cir- 
cularization radius; we adopt the narrower tidally trun- 
cated profile. For simplicity we assume a disk with uni- 
form temperature and surface gravity and hence no vari- 
ation in the emergent spectrum with location. This is a 
gross simplification but probably no worse an approxima- 
tion than the one we make by using stellar atmosphere 
spectra in place of more realistic disk atmospheres. The 
temperature distribution in quiescent disks is in any case 
expected to b e rather flatter than in outbursting systems 
l)Lasotal 120011 and references therein). Note that since 
we are calculating line profiles in the companion star's 
rest frame, the disk appears to move with velocity semi- 
amplitude K\ + K2 (where K\ and K% are the velocity 
semi-amplitudes of the compact object and companion 
star respectively). 

We show model composite line profiles in Fig. ^ as- 
suming equal fluxes from star and disk, and lines of the 
same equivalent width. Examining the high resolution 
case first, it can be seen that the much higher rotational 
broadening in the disk has the effect of completely sepa- 
rating the two components, even though the local emer- 
gent lines were assumed to be identical. The effect is 
amplified even more when the orbital average is consid- 
ered, as this further smooths the disk profile out. This 
separation is well preserved for a typical optical reso- 
lution (R ~ 4 000), and the two components are suffi- 
ciently decoupled that the apparent depth of the optical 
absorption line is little affected. The disk line is effec- 
tively broadened by so much that it appears almost as a 
smooth continuum on the scale of the stellar features. At 
even lower resolutions the separation becomes less clear, 
however. There are obviously still two components, but 
they merge more and one might expect a more signifi- 
cant contamination. Nonetheless, if one could isolate a 
single profile at high enough signal-to- noise, the stellar 
component could be separated and measured reasonably 
reliably. In practice this is not what is done of course, 
and the broad disk components would prove more diffi- 
cult to separate when multiple blended lines are present 
and when one must rely on statistical measures, such as 
X 2 fitting, to estimate the veiling. 

3.3. Synthetic Spectra 

We next convolve the above line profiles with model at- 
mosphere spectra. We construct these using atomic and 
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Fig. 2. — Simulated spectra of a K4V star, representative of the 
companion to A0620-00, and a disk with the same temperature 
and gravity. This illustrates how severely the spectrum is modified 
by the disk's rotational broadening (in the disk rest frame), and 
orbital smearing (in the stellar rest frame). Only the strongest 
features remain identifiable, and the distinctive double-troughed 
disk line profile is only visible in the disk rest frame. Successive 
vertical offsets of 0.2 units have been applied. 



molecular line data and a grid of local thermodynamic 
equilibrium (LTE) ATLAS9 model atmospheres obtained 
from the Kurucz website 3 . We interpolate within these 
models for a given set of T e g, logg, and [Fe/H] and then 
compute mod el spectra usin g the LTE stellar synthesis 
code, MOOG l|Snedenlll973T) . 

Fig.[3shows examples of stellar and disk spectral mod- 
els calculated in this way for the same assumed temper- 
ature and gravity, to illustrate the gross differences be- 
tween the two spectra, and reinforce the point made in 
the previous section. 

4. SYSTEMATIC ERRORS IN VEILING ESTIMATES 
4.1. Companion temperature 

A serious problem with the veiling estimates arises 
from the uncertainty in the spectral type. These esti- 
mates are extremely sensitive to a mismatch between the 
te mplate and the data, as is illust r ated f or GS 2000+25 
by Harlaf tis. Home, fc Filippenkol (^996), for example, 
where disk contributions of 0-35% are inferred for Kl- 
M0 templates. In theory, the best fitting template will 
indicate the true spectral type, but the correct choice 
may not be obvious from \ 2 values alone, as is the case 
for the GS 2000+25. 

To further explore this problem we perform a series of 
simulations in which we fit the T e g = 4690 K, log g = 4.46 
template plus featureless disk to composite models con- 
taining stellar spectra of a range of temperatures plus 
a featureless disk. In each case, we adjust the assumed 
(true) veiling, a, until the derived veiling, /3, matches 
that observed. Fig.[3]shows the assumed veiling required 
for a range of stellar temperatures roughly correspond- 
ing to the temperature uncertainty in A 0620-00. Large 
errors are clearly possible: if the true temperature is just 
100 K smaller than that of the template then the opti- 
cal veiling is underestimated by a factor of 2, and if it 

3 http://kurucz.harvard.edu 
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Fig. 3. — True veiling, a, required to reproduce observed values 
if the template is a poor match in temperature to the star. Cal- 
culations assumed an observed optical veiling of 6% (solid line), 
and an upper limit of 27% in the IR (dashed line), and a 4690 K 
template. 



is 100 K larger then no veiling is req uired at all. In- 
deed [GonzSe^Hemdnde^^T] l)2004|) derived zero red 
veiling in A 0620-00 for their preferred T cS = 4900 K 
solution. If the measured veiling is larger then the frac- 
tional error is less significant, but will still usually ex- 
ceed the statistical error, and will probably always dom- 
inate the uncertainty in veiling measurements. Similar 
e ffects occur in the IR, where the quoted up per limit 
of Shahbaz, Band vopadhvav. fc Char lea l)1999|) could be 
consistent with a true veiling of up to 50% for severe 
template mismatch. 

In both the optical and IR, if the temperature of the 
template is lower than that of the star, the veiling will 
be overestimated because cooler templates generally have 
stronger lines (averaged over the wavelength region being 
fitted). In the case of A 0620-00, there are indications 
that that may be the discussed above, so existing 

veiling estimates are probably useful as upper limits. In 
most systems, however, we do not have such a precise 
temperature and cannot make such an assumption. 

4.2. Companion gravity 

As noted above, the uncertainty in gravity for A 0620- 
00 is not large. Nonetheless, the CO bands in partic- 
ular are expected to be sensitive to gravity, so even a 
small uncertainty might be significant. Furthermore, the 
gravity of A 0620-00, and most BHXRT companions, is 
lower than that of a main sequence st ar, and hence lower 
than typical templates. For example iGravl l)1992f) gives 
loggt = 4.60 for a K4V star. 

As above, we take a T c g = 4690 K, logy = 4.46 tem- 
plate plus featureless disk model and use it to fit syn- 
thetic optical spectra comprising a star of the same tem- 
perature but different gravities, plus a 6 % featureless 
disk contribution. An uncertainty in log g of 0.04 has 
negligible effect, and allowing a range of ±0.1 still yields 
derived ve ilings of 6 ± 2 %, within the statis tical uncer- 
tainties of iMarsh. Robinson, fc Wood! l|1994D . The most 
extreme case is when we use a main-sequence log g — 4.60 
template and fit to a log g = 4.35 spectrum. Even then 



a derived veiling of 6 % corresponds to a true veiling of 
2%, close to the edge of the quoted error range. There- 
fore uncertainties in gravity, and mismatch with a main- 
sequence template, are not the dominant source of error 
in the optical. 

In spite of the gravity sensitivity of the CO bandheads, 
there is also negligible effect in the IR. Using a logg = 
4.46 template, we a mismatch in gravity only raises the 
upper limit from 27% to 29% for gravities of log g = 
4.46 ± 0.10. Fitting a \ogg = 4.60 template to a model 
spectrum with gravity of log g = 4.35 only raises the limit 
to 30 %. Note that the bias is in the opposite sense to 
that found in the optical. 

4.3. Companion metallicity 

Another difficulty is that the line strengths will also 
be a function of the metallicity. A mismatch between 
the metallicities of the template and companion could 
also introduce an error in the veiling measurement; a low 
metallicity companion will have weaker lines, and hence 
appear more heavily veiled than it actually is. 

The metallicity of the secondaries of BHXRTs is 
not w ell constrained by observations. iKotoneva et alJ 
(2002) examine metallicities of nearby K dwarfs and 
find a spread in [Fe/H] of approximately —1.0 to 
+0.2 . The distribution peaks at -0 .2. For compari- 
son, IGonzalez Hernandez et all l)2004|) derive [Fe/H] = 
+0.14 ± 0.20 for A 0620-00, at the high end of the 
range. Given this range, unless the template star 
is deliberately selected by metallicity, a mismatch is 
not just possible, but likely . The template used 
by IMarsh. Robinson, fc Wood! lll99l. HP 16160, has 
[Fe/H] = -0.03 + 0.14 l|Heiter Lucid 12001 . so 
the mismatch is small, but not negligible (0.17 ± 
0.24). In the case of the IR measu rement of 
Shah baz. Bandvopadhvav. fc Charlesl (^999) , however, 
the preferred template (HD 42606) does not have a pub- 
lished metallicity, so a larger mismatch is possible. 

To test the sensitivity to metallicity we adopt a model 
spectrum of T e s = 4690 K, log g = 4.46, solar metallic- 
ity, and a variable disk contribution. The template we fit 
to the model is assumed to have the same temperature 
and gravity, but metallicities o f -0.3, -0.6, and - 0.9, i n 
accord with the range found bv IKotoneva et all |2002). 
With these combinations we find that to reproduce an 
observed optical veiling of 6 % requires true veilings of 
20 %, 32 %, and 46 %. In the IR the sensitivity to metal- 
licity is also large; an upper limit on the veiling of 27% 
permits true veilings of up to 41%, 54%, and 67% for 
template metallicities of -0.3, -0.6, and -0.9 respectively. 
Thus if the template has lower metallicity than the tar- 
get, then the veiling could have been substantially un- 
derestimated. 

Ret urning to A0620-00, IGonzalez Hernandez et all 
1)20041) obtain a metallicity of \Fe/H] = 0.14 ± 0.20 
where as the template used bv lMarsh. Robinson, fc Wood! 
(1994) has a metallicity of [Fe/H] = -0.03 ± 0.14. On 
the one hand the difference between them is not sta- 
tistically significant, but on the other it could be large 
enough to bias the veiling measurement substantially, by 
a factor of two or more. Mismatch between template 
and target metallicities can thus be a major contributor 
to the uncertainty in veiling measurements. 
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Fig. 4. — Effect on veiling estimates if the disk has absorption 
lines. The upper panel uses a red spectrum (5920-6520 A), and 
the lower panel a section of the K band (2.19—2.31 )im). This 
is for an ideal case where the template spectrum is the same as 
was used in the simulation. Even in this case, significant errors 
occur, especially at low resolution, or when broad bands such as 
CO dominate as is the case in the K band. 



4.4. Disk spectral features 

We finally proceed to test the robustness of veiling esti- 
mates if the disk spectrum is not featureless. We consider 
an ideal case, where the spectra are averaged over a full 
binary orbit in the rest frame of the companion. We use 
the modeled companion star spectrum (T c s — 4690 K, 
\ogg = 4.46) as the template for the stellar component 
so that no bias is introduced by template mismatch. Val- 
ues of /3/a are shown in Fig. 0] for various assumed disk 
temperatures and resolutions; we have used a = 0.1 for 
these tests. 

We consider first the optical case. For moderately high 
resolutions (R — 4, 000 and R — 30, 000), the error in the 
veiling is modest, and would usually be within the sta- 
tistical uncertainties in the veiling measurement (which 
are typically large). The error is larger for a cool disk 
because the spectral features tend to be deeper. At low 
resolutions, the problem becomes much more severe. In 
this case, as can be seen in Fig.^ the disk and stellar line 
profiles are no longer well resolved, and it becomes hard 
for the fitting algorithm to distinguish between them, 
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even with a template perfectly matched to the star. For 
low disk temperatures, the stronger lines actually mean 
that adding a disk contribution makes the overall lines 
appear stronger rather than weaker, and one would de- 
rive a negative veiling (clipped here to zero). 

In the IR the situation is less favorable, at least if the 
primary diagnostics are CO bands. Data are typically 
obtained at lower resolution than in the o ptical (as in 
Shah baz. Bandvopadhvav. fe Ch arles 1999}), and even if 
higher resolutions are used, the intrinsic width of the CO 
bands limits the discrimination possible between disk and 
stellar profiles. The low resolution case is thus compara- 
ble to the same resolution in the optical, but there is less 
benefit to high resolutions, unless the study focuses only 
on metallic features. As noted earlier, contamination of 
the CO bands by the disk might be expected in quiescent 
BHXRTs, as CO emission is seen in the quiescent cat- 
aclys mic variable WZ Sge l|Howell. Harrison, fc Szkodvl 
I2004J) . although it is not clear that the cases are compa- 
rable. If contamination is by emission features then the 
veiling will tend to be overestimated. 

We thus conclude that provided veiling studies are 
done at intermediate resolution or better, the bias on 
the veiling is not the dominant source of error in the op- 
tical, even if the disk has strong absorption lines. At low 
resolutions, such as commonly used in the IR, the mea- 
surement could be severely compromised. This will also 
apply (irrespective of resolution) if the features used are 
not narrow, e.g. CO or TiO molecular bands. 

5. EFFECT ON THE BROAD-BAND SPECTRAL ENERGY 
DISTRIBUTION 

It is usually argued (e.g. iGelino. Harrison, fc Oroszl 
12001 that since the veiling is measured to be lower in the 
red than in the blue, then the IR veiling must be even 
lower. The implicit assumption being made here, how- 
ever, is that the veiling varies monotonically with wave- 
length. There is no requirement that this be the case; 
multiple components may contribute to the disk spectral 
energy distribution. The optical veiling (and UV spec- 
trum) may well be dominated by a stream impact point 
or overflow l|McClintock et al.ll2003j) . and appears to be 
relatively blue. The stream-impact point certainly does 
appear as a bright source of emission line light in Doppler 
tomograms l)Marsh. Robinson, fc Woodl I 19941k In con- 
trast, it is observed that rapid flaring exhibits a rather 
redder sp ectrum, comparable to the color of the compan- 
ion star (Sha hbaz et alJl2008l) : this could represent the 
optically thin disk component. Optically thick disk ma- 
terial could also be present. This material is predicted 
to be very cool (~3000K) and will not contribute much 
flux at optical wavelengths. Given the large projected 
surface area of the disk (several times that of the com- 
panion), an extended but cool component might, how- 
ever, produce a significant contribution in the IR, even if 
it is undetectable in the optical. One can then imagine a 
disk with two components, a hot component dominating 
optical and UV wavelengths and a cool component domi- 
nating infrared wavelengths. Such a picture is consistent 
with theoretical expectations, but has the unhappy con- 
sequence that the total veiling may well be minimized in 
the red, precisely where most measurements are made! 

We can attempt to quantify the possible contribution 
of the cool disk material to the IR by considering opti- 
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Fig. 5. — Upper limits on the temperature and covering frac- 
tion (relative to a disk extending to the tidal truncation ra- 
dius) of an optically thick disk in A0620-00. The solid line 
indicates constraints due to th e observing veiling in the red 
IGonzal cz Hernandez ct al. 2004 t). Dotted lines are base d on 
the veiling at K iShahbaz, Bandvopadhvav. & Charles 1999) and 
dashed lines are based on t he consistency of inclination determina- 
tions across J, H, and K IGelino. Harrison, fc O rosz 2001}) . The 
veiling limits assume a K4V spectral type as this places stricter 
limits than earlier types. Conversely, the inclination constraint is 
based on K2V, as this is stricter than later types. Veiling con- 
straints correspond to 2cr limits (~ 90% confidence). The inclina- 
tion constraint uses 1 a limits, but is effectively based on a joint 
confidence region in the J and K (with their errors in opposite 
directions), so is also close to 90% confidence. 



cally thick disks with a range of temperatures and pro- 
jected areas. As above, we assume a uniform disk tem- 
perature and parameterize the area as the covering factor 
of a tidally truncated disk. For A0620-00 such a maximal 
disk has a projected area ~ 7 times that of the compan- 
ion st ar, assuming i — 40? 75 l)Gelino. Harrison. fc Oros3 
|2001[K whereas a disk only extending to the circulariza- 
tion radius still subtends an area ~ 3 times that of the 
companion. We can then estimate the maximum cov- 
ering factor which is consistent with the observed opti- 
cal veiling as a function of temperature. For this pur- 
pose we now use the most robust estimate , a re d veil- 
ing of ± 5% (|Gonzalez Hernandez et aLll200l . We 
estimate the relative contributions from the companion 
and disk as a function of temperature by using surface 
brightness from NextGen model atm osphere calculations 
l)Hauschildt. Allard. fe Baronl ll999). This of course ex- 
plicitly assumes the disk is optically thick, so these con- 
straints will not apply if the disk is optically thin. 

We find that disks of a temperature comparable or 
greater than the companion are highly constrained by 
the optical veiling if this is reliable (see Fig.[SJ); this con- 
straint would become even tighter if one allows for the 
presence of a blue veiling component as well. What this 
means is that the quiescent disk must either be much 
cooler than the companion (< 3500 K) or mostly opti- 
cally thin (with perhaps small clumps of optically thick 
material). Note that since the circularization radius ef- 
fectively limits the minimum size of the disk, one cannot 
explain this away entirely with a small disk; the mini- 
mum covering factor possible in this way is about 0.4. 

As the disk temperature is allowed to drop, how- 



ever, the constraint relaxes, and for theoretically m o- 
tivated temperatures (~3 000K; e.g. iMenoul [20021) a 
more respectable covering fraction is possible. Such a 
disk can actually contribute strongly in the K band 
while having negligible effect in the optical. It would 
of course veil the features of the companion star, so 
K band veiling estimates could still constrain its pres- 
ence, but as we have seen the existing I R estimate 
(Shahbaz. Bandvopadh vav. fc Charles! 1199 9^1 is poten- 
tially compromised by even a small mismatch with the 
template temperature or metallicity, or by disk spectral 
features. The upper limit derived, 27 % at 2-cr, however, 
does not strongly constrain a cool disk (Fig. 03), especially 
when we consider that the upper limit may not be very 
robust. 

An additional constraint comes in the form 
of the consistency of inclination determinations 
from the simultaneous J, H , and K lightcurves 
l)Gelino. Harrison, fc Oroszll200H Gelino 2005, private 
communication). These authors find inclinations de- 
rived from individual bands are all consistent with the 
combined fits to within uncertainties (±1.5°). If veiling 
were present with colors different to those of the star, 
then different bands would experience different veilings 
and hence discrepant inclinations would be derived. We 
show in Fig. 03 the constraint this places on the disk 
covering factor with a K2V companion star. Assuming 
K4V results in slightly higher limits. Unfortunately 
even these observations only limit the covering factor, 
and hence veiling, to a similar precision as claimed from 
spectroscopic veiling. 

A final possible constraint would come from the ob- 
served spectral energy distri bution (SEP). Indeed one 
of the arguments made by IGelino. Harrison, fc Oros3 
l)2001[) for negligible IR contamination was that the pho- 
tometric SED was consistent with the spectroscopic spec- 
tral type constraints. This is less robust than the inclina- 
tion constraint, however, as it is very sensitive to the as- 
sumed spectral classification of the companion star, and 
hence results in weaker limits than the inclination con- 
straint does. These have not been shown. 

We conclude that an extensive, cool, optically thick 
disk component in A0620-00 cannot confidently be re- 
jected based on any existing observational evidence. 
There is no reason such a disk should not be present, 
and it could have as large a contribution in the K band 
as veiling measurements allow, 25%, or even more. We 
know of no evidence to support the assertion that the K 
band contamination is as low as a few percent. 

6. DISCUSSION 

A cool, optically thick disk component might be ex- 
pected in quiescent BHXRTs as it has long been argued 
that an entirely optically th in disk is incon siste nt with 
the D IM (see discussions by lCannizzolll993l and iMenoii 
l20f)l . The argument is made that an optically thin 
disk requires extremely l o w sur face densities, for example 
iLin. Williams, fc Stove! l|1988fl estimated S < 1 gem" 2 . 
In contrast, the DIM oper ates at critical surface densi- 
ties o f order 100 gem -2 (e.g. lHameurv et alll998llLasotal 
2001). Even if one beli eves that the DIM is not cor- 
rect. iMenou et aTl l)2000l) estimate that during a typical 
outburst of A0620-00 a few x 10 24 g are accreted, corre- 
sponding to an average surface density pre-outburst of 
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at least 20 gem . 

While it would appear that most of the quiescent disk 
mass must then be in the form o f optically thick material, 
there is an inconsistency here. iLin. Williams, fc Stoverl 
(1988) were considering relatively high temperature ma- 
terial (by quiescent standards) at temperatures above 
6000 K, whereas disk models are now working with tem- 
peratures that drop to 2000-3000 K in quiescence. At 
these temperatures opacit ies are expected to be lower 
(e.g. iFerguson et aLlfeOOSj) . so higher surface densities 
could remain optically thin under these conditions, and 
this issue requires further examination. 

As a crude test we considered the expected optical 
depth, based on Ros seland mean opacities taken from 
IFerguson et alJ l)2005j) of a range of uniform slabs of tem- 
perature 2000-3000K, surface density 20-100gcm- 2 , 
and thickness l-4xl0 9 cm, corresponding approximately 
to H/R values of 0.01-0.02 for the outer disk of a typical 
quiescent BHXRT such as A 0620-00. At low tempera- 
tures (T ~ 2000 K) this essentially depends only on sur- 
face density, with optical depths 0.2-1.2 likely. At higher 
temperatures (T ~ 3000 K) volume density becomes im- 
portant, and a larger range of optical depths from 0.1- 
3.8 are found for the parameters described above. We 
conclude that based on current knowledge of such disks 
optically thick regions might or might not be present; 
they do not appear to be inescapable as has sometimes 
been claimed, but should still be considered a possibility, 
and observational constraints such as discussed here are 
needed to rule it out. 

As has long been appreciated, an optically thick disk is 
not inconsistent with the presence of disk emission lines, 
as these likely originate from a chromosphere above the 
disk; indeed this chromosphere could account for most or 
all of the optical emission. Alternatively it may be that 
there are optically thin p atches in the disk in between op- 
ticall y thick clumps (e.g. iVrielmann. Hessman. fc Hornd 
2002). Such a picture would be appealing to explain 
the mirror eclipses seen in quiescent dwarf novae, which 
imply that at least 20 % of the disk (by radius) mus t be 
opti cally thin in the cont inuum l|Littlefair et all2001l: see 
also lFroning et al]ll999jh Note that if the disk mass is 
clumpy then this will further increase the likelihood of 
clumps being optically thick. 

Any optically thick material must be cool. For it to be 
stable on the cool branch of the disk instability requires 
T e g < 6000 K. Unless it consists of a few very dispersed 
clumps, or is indeed optically thin, the optical veiling 
further constrains the temperature to T c g < 3500 K. 
This is consistent with typical t heoretical temperatures 
T cB ~ 3000 K fe.g. lMenoiJ l20021. In this regime, as we 
have seen, the observational constraints are extremely 
weak, and an extensive optically thick component con- 
tributing a significant fraction of the light in the K band 
is possible. 

It would therefore be an interesting constraint if the 
disk contribution at K is only a few percent. As shown 
in Fig. El this would require the area of the optically 
thick fraction of the disk to be small, implying most of 
the disk is optically thin. A significant optically thick 
covering factor (more than 10 %) is only consistent with 
5% contamination for T e g < 1800 K. It is unlikely that 
quiescent disk material can be this cool, as there will al- 
ways be some heating by the companion star, the stream 



impact, and the central X-ray source. 

One argument against significant IR contamination is 
that optical ellipsoidal modulations are distorted and 
variable, but the IR modulations are much cleaner. This 
argument is not conclusive, however. It is likely that the 
optical veiling originates from the stream impact and/or 
overflow. This material is highly non-axisymmetric, and 
variable. The optically thick material, however, is likely 
to be in an axisymmetric and probably non- variable disk, 
and so is likely to only contribute a steady source of light 
that does not distort the ellipsoidal variations. 

7. CONCLUSIONS 

We have examined the observational constraints upon 
the presence of a cool, optically thick component to the 
accretion disks of quiescent BHXRTs. Such a component 
is theoretically plausible, but has not yet been detected. 
We have shown that constraints from veiling of the stel- 
lar photospheric spectrum are not very robust, as the 
veiling derived is very sensitive to a mismatch between 
temperatures and/or metallicities of the companion star 
and template, and to a lesser extent to a surface gravity 
mismatch. If the disk spectrum exhibits similar spectral 
features to the star, for example CO bands in emission or 
absorption, then this can further bias veiling estimates. 

The current situation in the prototypical system 
A 0620-00 is that reasonably robust optical veiling mea- 
surements now exist, in the sense that the veiling, 
temperature, gravity, and metallicity are all deter- 
mined simultaneously by high r esolution spectroscopy 
l)Oonzalez Hernandez et all 12004^ . but that the IR re- 
mains poorly constrained. The red veiling tightly con- 
strains the temperature and or covering factor of a warm 
disk, requiring either T c ff < 3500 K, or a very low cover- 
ing factor (< 10 %) of optically thick material. Cooler 
disks, as expected theoretically, are unconstrained by 
the optical measurement. The IR veiling measurements 
made so far are sensitive to the uncertainties summarized 
above so are not robust. A more robust approach yielding 
similar constraints to the IR veiling is the relative ampli- 
tudes of ellipsoidal variations l)Celino. Harrison, fc Oroszl 
120011) . but these still allow contamination up to ~ 25%. 

Thus we conclude that no current measurements con- 
strain the K band contribution from the accretion disk 
to better than about 25 % in A 0620-00, and conse- 
quently that significant dilution of IR ellipsoidal mod- 
ulations remains possible. In particular, the assump- 
tion that the IR contamination is only a few percent is 
not yet supported by compelling evidence. This issue 
has imp ortant implications for black hole mass determi- 
nations. iGelino. Harrison, fc Oros3 l)2001fl . for example 
commented that increasing the K band disk veiling to 
50% would increase the derived inclination to 60°. The 
derived black hole mass would then drop from 1 1 M Q to 
below 5M©! This is an extreme case, clearly a milder 
and more plausible IR veiling could still have substantial 
impact on the derived black hole mass. 

Finally we note that if the K band disk contamina- 
tion can be observationally constrained to be only a few 
percent then this will place very tight constraints indeed 
upon the temperature and/or covering factor of optically 
thick material, providing a test of quiescent disk models. 

RIH is grateful to Dawn Gelino for helpful criti- 



Disks in Quiescent Black Hole Binaries 



9 



cism of this manuscript, to Rick Hessman for discus- 
sions of patchy accretion disks in quiescent objects, 
to Tariq Shahbaz on optimal subtraction techniques. 
RIH acknowledges support for most of this work from 
NASA through Hubble Fellowship grant #HF-01150.01- 



A awarded by the Space Telescope Science Institute 
(STScI), which is operated by the Association of Univer- 
sities for Research in Astronomy, Inc., for NASA, under 
contract NAS 5-26555. 



REFERENCES 



Al-Naimiy, H. M. 1978, Ap&SS, 53, 181 

Cannizzo, J. K. 1993, in Accretion Disks in Compact Stellar 

Systems, Ed. J. C. Wheeler, 6 
Cherepashchuk, A. M. 2000, Space Science Reviews, 93, 473 
Ergma, E., & Sarna, M. J. 2001, A&A, 374, 195 
Ferguson, J. W., Alexander, D. R., Allard, F., Barman, T., 

Bodnarik, J. C, Hauschildt, P. H., Heffner-Wong, A., & 

Tamanai, A. 2005, ApJ, 623, 585 
Froning, C. S., Robinson, E. L., Welsh, W. F., & Wood, J. H. 1999, 

ApJ, 523, 399 
Froning, C. S. & Robinson, E. L. 2001, AJ, 121, 2212 
Gelino, D. M., Harrison, T. E., & Orosz, J. A. 2001, AJ, 122, 2668 
Gonzalez Hernandez, J. I., Rebolo, R., Israclian, G, Casares, J., 

Maeder, A., & Meynet, G. 2004, ApJ, 609, 988 
Gray, D. F. 1992, Cambridge Astrophysics Series, Cambridge: 

Cambridge University Press, 1992, 2nd ed. 
Gray, R. O. & Corbally, C. J. 1994, AJ, 107, 742 
Hamcury, J., Menou, K., Dubus, G, Lasota, J., & Hure, J. 1998, 

MNRAS, 298, 1048 
Harlaftis, E. T., Home, K., & Filippenko, A. V. 1996, PASP, 108, 

762 

Harrison, T. E., Osborne, H. L., & Howell, S. B. 2004, AJ, 127, 
3493 

Haswell, C. A. 1992, Ph.D. Thesis, The University of Texas at 
Austin 

Haswell, C. A., Hynes, R. I., King, A. R., k, Schenkcr, K. 2002, 

MNRAS, 332, 928 
Hauschildt, P. H., Allard, F., & Baron, E. 1999, ApJ, 512, 377 
Heitcr, U. & Luck, R. E. 2003, AJ, 126, 2015 

Howell, S. B., Harrison, T. E., & Szkody, P. 2004, ApJ, 602, L49 
Hynes, R. I., Zurita, C, Haswell, C. A., Casares, J., Charles, P. A., 

Pavlenko, E. P., Shugarov, S. Y., & Lott, D. A. 2002, MNRAS, 

330, 1009 



Hynes, R. I., Charles, P. A., Casares, J., Haswell, C. A., Zurita, C, 

& Shahbaz, T. 2003, MNRAS, 340, 447 
Hynes, R. I., et al. 2004, ApJ, 611, L125 

Kotoneva, E., Flynn, C, Chiappini, C, & Matteucci, F. 2002, 

MNRAS, 336, 879 
Lasota, J. -P. 2001, New Astronomy Review, 45, 449 
Lin, D. N. C, Williams, R. E., & Stover, R. J. 1988, ApJ, 327, 234 
Littlefair, S. P., Dhillon, V. S., Marsh, T. R., & Harlaftis, E. T. 

2001, MNRAS, 327, 475 
McClintock, J. E. & Remillard, R. A. 1986, ApJ, 308, 110 
McClintock, J. E., Narayan, R., Garcia, M. R., Orosz, J. A., 

Remillard, R. A., & Murray, S. S. 2003, ApJ, 593, 435 
Marsh, T. R., Robinson, E. L., & Wood, J. H. 1994, MNRAS, 266, 

137 

Menou, K. 2002, ASP Conf. Ser. 261: The Physics of Cataclysmic 

Variables and Related Objects, 387 
Menou, K., Hameury, J., Lasota, J., & Narayan, R. 2000, MNRAS, 

314, 498 
Oke, J. B. 1977, ApJ, 217, 181 

Shahbaz, T., Bandyopadhyay, R. M., k, Charles, P. A. 1999, A&A, 
346, 82 

Shahbaz, T., Dhillon, V. S., Marsh, T. R., Zurita, C, Haswell, 
C. A., Charles, P. A., Hynes, R. I., & Casares, J. 2003, MNRAS, 
346, 1116 

Sneden, C. 1973, ApJ, 184, 839 

Vrielmann, S., Hessman, F. V., & Home, K. 2002, MNRAS, 332, 
176 

Zurita, C, Casares, J., & Shahbaz, T. 2003, ApJ, 582, 369 



